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LS 5039 is a gamma-ray binary system observed in a broad energy range, from radio to TeV 
energies. The binary system exhibits both flux and spectral modulation as a function of its orbital 
period. The X-ray and very-high-energy (VHE, E > 100 GeV) gamma-ray fluxes display a maxi¬ 
mum/minimum at infer!or/superior conjunction, with spectra becoming respectively harder/softer, 
a behaviour that is completely reversed in the high-energy domain (HE, 0.1 < E < 100 GeV). The 
HE spectrum cuts off at a few GeV, with a new hard component emerging at E > 10 GeV that is 
compatible with the low-energy tail of the TeV emission. The low 10 - 100 GeV flux, however, 
makes the HE and VHE components difficult to reconcile with a scenario including emission 
from only a single particle population. We report on new observations of LS 5039 conducted with 
the High Energy Stereoscopic System (H.E.S.S.) telescopes from 2006 to 2015. This new data 
set enables for an unprecedentedly-deep phase-folded coverage of the source at TeV energies, as 
well as an extension of the VHE spectral range down to ~120 GeV, which makes LS 5039 the 
first gamma-ray binary system in which a spectral overlap between satellite and ground-based 
gamma-ray observatories is obtained. 
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1. Introduction 

LS 5039 is a binary system located at ~ 3.5 kpc from the Earth, only visible from the Southern 
hemisphere. It is composed of a compact object in orbit around a 06.5V star, with an orbital 
periodicity of Porb = 3.90603 days [1], a moderate eccentricity (e = 0.35), and an angle position in 
the plane of the sky barely constrained (/« 20 — 65°). The mass of the companion star is ~ 23 Mq, 
and its radius 9.3 Rq (Fig. 1). The nature of the compact object (Mx ~ 3.7 Mq), either a black hole 
or a neutron star, remains unclear. 

LS 5039 has been detected in radio, featuring persistent outflows (with sizes ranging ~ 2 — 
2000 astronomical units, AU) which classified it as a possible microquasar candidate [2]. In X-rays 
and high-energy (HE) gamma-rays, the emission appears modulated with the orbital period [5]. At 
very high energies (VHE), the High Energy Stereoscopic System (H.E.S.S.) telescopes detected 
the system in 2004 with ~ 11 h of observations obtained during the H.E.S.S. Galactic Plane survey 
[3]. After further follow-up on the source, a modulation of its VHE gamma-ray flux was revealed, 
allowing for a determination of the periodicity in gamma-rays, with P = 3.9078 ± 0.0015 d [4]. 

Markedly different flux and spectral properties were observed close to the system’s superior 
and inferior conjunction (SUPC and INFC, respectively). A lower gamma-ray flux and a spectrum 
well-fit by a simple power law (PL) is derived for SUPC. At INFC, the source is brighter and the 
spectrum significantly deviates from a power law above a few TeV Both hardness and intensity 
at TeVs seem anti-correlated with what is obtained at HEs [5]. Such anti-correlation can be ef¬ 
fectively explained by the angular dependence of both inverse Compton (IC) and pair-production 
cross sections and the variable observing conditions of the system along the orbit. However, the 
properties and location of both accelerator(s) and emitter(s) in the source are far from understood 
(see, e.g., [6, 7, 8, 9, 10, 11], and references therein). 

Since the last H.E.S.S. publication in 2006 [4], a substantial amount of observation time has 
been devoted to LS 5039. In the next section we report on an updated analysis of H.E.S.S. phase 
I data, whereas Sect. 3 focuses on the first results obtained with the H.E.S.S. phase II, for which 
the so-called mono analysis [15] (the new 28 m telescope, CT5 hereafter, in standalone mode) has 
been employed. A summary and main conclusions of this study are presented in Sect. 4. All results 
shown in this proceedings, albeit thoroughly cross-checked with an independent data calibration 
and analysis framework, are to be considered preliminary. A more detailed report is in preparation. 


2. Observations with H.E.S.S. phase I 

New data on LS 5039 have been obtained from 2006 to 2012 with the H.E.S.S. phase I tele¬ 
scopes. In parallel, more sophisticated analysis techniques and H.E.S.S. data-reduction software 
have been developed, notably improving both the background rejection power, the accuracy on 
the energy reconstruction of the gamma-ray showers and the angular resolution of the instrument 
[12, 13]. Below we report on the results of this new data set, as well as a re-analysis of pre-2006 
observations with the latest analysis tools, to characterise the long-term behaviour of LS 5039 at 
VHEs in a time-range spanning more than 8 years. The updated H.E.S.S. I data set amounts to 
about ~ 104 h of live observation time, with a mean zenith angle of 19.8°, yielding a detection of 
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Figure 1: Geometry of the orbit (in the orbital plane) of the compact object around the 06.5V star (Casares 
et al. [1]). The star size is to scale. The arrow at the bottom indicates the direction to the observer. Sev¬ 
eral notable positions of the compact object are indicated: periastron (0 = 0), apastron (0 = 0.5), superior 
conjunction (0 = 0.058) and inferior conjunction {(p = 0.716). 


the source at a statistical significance of more than 56a and 2800 detected gamma-rays. To perform 
this analysis, we used a ring background subtraction [16] in stereo mode. 

2.1 H.E.S.S. I: phase-folded light curve 

The integral flux above 1 TeV, folded with the orbital period of the system, is computed as¬ 
suming a spectral index F = 2.20 ± 0.03, where F is derived from a fit with a simple power law, 
dN/dEy oc Ey^, to the whole data-set. In Fig. 2 the known modulation of the source at VHE is 
recovered, which keeps essentially unchanged from one year to the next, for the whole data set, 
averaging over short, e.g., daily time-scales. 

To better characterise this modulation. Fig. 3 displays the phase-folded averaged flux for phase 
intervals of width 0.1 each. In this case, however, the large data set obtained in the last ~ 10 yr 
of observations enables us to fit separately each phase interval to compute the bin-averaged flux. 
We assume here again a simple power law model fitted to the data in each phase interval. The 
corresponding photon index F plotted against the differential gamma-ray flux at 1 TeV in the same 
bin is displayed in Fig. 4. We notice a relationship between both variables, the correlation coef¬ 
ficient is —0.93. When the system is brighter the spectrum is harder. In the same way, when the 
differential flux decreases, the spectral index increases. We note that noticeable departures from a 
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Phase 

Figure 2: Run-by-mn light curve corresponding to pre-2006 data (red), the same data set as reported in 
[4], reanalysed here with the same analysis tools as for the post-2006 data set (blue). The horizontal upper, 
middle and lower dotted lines represent the average flux for the INFC, all-data, and SUPC, respectively. 



Figure 3: Phase folded light curve, with LS 5039 integral flux above 1 TeV plotted over the system orbital 
phase. Integral fluxes are computed from a separate power-law fit in every phase interval. 

power law may not be excluded for some phase intervals. A more detailed spectral characterisation 
is in preparation and will be shown in a future publication. 

The light curve in Fig. 3 displays a ratio > 8 between the maximum and minimum flux levels, 
corresponding to the compact object being at the INFC and SUPC, respectively, with an absolute 
maximum at (^ ~ 0.75 ± 0.05. A detection of LS 5039 at a significance level above 5a is now 
obtained in every single phase-bin. We did not try, however, to fit the light curve neither to constrain 
its timing structure nor to characterise its asymmetry or the presence of one or more components, 
which is left for a forthcoming publication. 

2.2 H.E.S.S. I: spectral analysis 

The large data set now available allows for a rather high-precision characterisation of the 
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Figure 4: Differential gamma-ray flux at 1 TeV versus photon index F obtained from a dedicated power-law 
fit in each 0.1-width phase interval. The dotted line and the statistics in the legend correspond to a linear fit 
to the data. 



Figure 5: Spectral energy distribution for the analysis of H.E.S.S. 1 data obtained in 2004 - 2012 The data 
have been divided between orbital phases close to SUPC (blue) and INFC (red), which are well fit by a 
simple power law and an exponential cutoff model, respectively. 


averaged spectra at both INFC (at orbital phases (j) G [0.45-0.9]) and SUPC (with (j) < 0.45 or 
(j) > 0.9), see Fig 5. A straight power law fits well the spectrum at SUPC, with an averaged spectral 
index of F = 2.406 ± 0.052. For orbital phases close to INFC, a power law with an exponential 
cutoff model dN/dEy oc x exp(—E^/Ecut) has been fitted (with F = 1.843 ± 0.063, and an 
energy cut Ecut = 6.6 ±1.6 TeV) in agreement with previous results [4]. Further spectral models, 
e.g. a broken-power law or a log-parabolic distribution are not reported here, although they also 
provide a better fit to the INFC data compared to a simple power law assumption. 
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Figure 6: LS 5039 excess map from a ring background method, for a mono analysis with ^ 19 h of ob¬ 
servation time and safe cuts, is displayed in the left panel. LS 5039 is clearly detected, together with the 
nearby, extended pulsar wind nebula HESS J1825-139. The distribution of the background events, after 
excluding both LS 5039 and HESS J1825-139, is displayed on the right panel, featuring a well-centered and 
well-normalised significance distribution. 

3. Observations with H.E.S.S. phase II 

3.1 H.E.S.S. 11: dataset and analysis used 

LS 5039 has been observed during the commissioning of H.E.S.S. II in 2012-2013, with further 
observations taken in 2014 and 2015. Here we present the first results of these observations, which 
amount to a total of 39 runs taken towards LS 5039 and yielding a total observing time of 18.5 
h of data after quality selection cuts. The data have been analysed in the so-called mono mode, 
making use of the newly deployed CT5 telescope in standalone mode, which permits to reach the 
lowest energy threshold attainable with H.E.S.S. II. An analysis performed within the so-called 
combined mode, with a similar performance at the lowest energies but with additional stereo-moAe. 
sensitivity above tens of TeV [15], will be presented in a forthcoming publication. This analysis 
makes use of a technique based on a semi-analytical shower development model [12], adapted to 
the CT5 mono mode [15] observations. All results have been cross-checked with an independent 
calibration software and analysis chain (ImPACT, [13], also adapted to CT5 mono observations), 
providing compatible results. 

The results reported here account for observations taken under zenith angles ranging in be¬ 
tween 9° and 52°, with an average of 31°. Data were taken in wobble mode, where the telescope 
alternates pointing directions with a given offset from the source position. We make use of the 
Ring background method [16] to produce sky-maps and significance distributions in the EoV (see 
left panel in Eig. 6. LS 5039 is clearly detected with at a statistical significance of 12a (with 
342 gamma-rays). In addition, the nearby source HESS J1825-139 is also clearly detected in this 
analysis, demonstrating the capabilities of this CT5 mono mode to reconstruct the signal from 
extended sources. The background significance distribution after subtracting both LS 5039 and 
HESS J1825-139 is well normalised, as displayed in the right panel of Eig. 6. 
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Figure 7: SED of LS 5039 including the analysis of H.E.S.S. phase-11 data in CT5 mono mode (black), 
observations taken between 2004 and 2012 with the H.E.S.S. phase-I telescopes (blue), and Fermi-LAOC 
spectral points (purple) extracted from Takata et al. (2014) [14], with a purple solid line representing the fit 
with a power law with an exponential cutoff in the energy range 0.1-10 GeV (with E = 2.06 ± 0.02 and 
Ecut = 3.42 ± 0.17 GeV). 


3.2 H.E.S.S. 11: spectral analysis 

We performed a CT5 mono spectral analysis in order to reconstruct events in the ~ 0.1 — 2 
TeV energy range. The Reflected background method, in which OFF events are taken from several 
regions at the same radial distance with respect to the camera centre as the target position, is 
employed to retrieve LS 5039’s spectrum. Standard cuts for the background subtraction have 
been used, allowing for a robust characterisation of LS 5039 spectra down to a threshold energy of 
119 GeV for the current data set analysis. Although lower energies may be reached after accounting 
for a larger H.E.S.S. data set, for which a dedicated study is in preparation, we note that this 
preliminary analysis already allows for the study of the energy range covered by both H.E.S.S. and 
the Fermi-LKI. This is the first time that a full GeV-TeV overlap is obtained for any gamma-ray 
binary system. 

A spectral energy distribution (SED) of ES 5039 at gamma-rays is displayed in Eig. 7. The 
spectral fit assumes a pure power-law model, from which a spectral index of F = 2.20 ± 0.03 
is retrieved, in perfect agreement with the analysis of (non contemporaneous) H.E.S.S. phase I 
analysis in the overlapping energy range. The H.E.S.S. phase II CT5 mono spectrum is compatible 
within the relatively large error bars of the > 100 GeV flux obfained in the analysis of Femii-LKF 
data reported in [14]. A refined analysis of fhe H.E.S.S. data, accounting for further observations 
of ES 5039 with CT5, as well as an updated Femii-LKF data set, will provide a deeper overlap of 
both instruments in a wider energy range, allowing for a more detailed study of the source spectral 
properties. 
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4. Summary 

The results reported here for the analysis of H.E.S.S. phase I data, both aeeounting for a 
reanalysis of previous observations with more sensitive analysis tools and the inelusion of new data 
taken from 2006 to 2012, confirms the regular behaviour of LS 5039 in a time range spanning more 
than ~ 8 yr. The evolution of both the VHE gamma-ray flux level and spectral properties depending 
on the orbital position of the compact object in the system are confirmed with this new analysis. 
A more detailed study, including phase-folded spectral characterisation in relatively thin phase 
ranges, the characterisation of the degree of contamination of the nearby source HESS J1825-139, 
and the implications in terms of emission/absorption processes responsible for the TeV emission, 
is in preparation. 

The analysis of the new H.E.S.S. phase II observations reported here allow for a spectral 
derivation down to ~ 120 GeVs for the data-set analysis results shown here. This is the first 
case for any known gamma-ray binary system in which a spectral overlap between satellite and 
ground-based observatories is obtained. Given the growing evidence that GeV and TeV emis¬ 
sion are produced in different regions and/or by separate particle populations, a deeper study of 
this overlapping energy range may provide strong constrains to the physics behind the gamma-ray 
emission in ES 5039 and possibly other similar systems. 

ACKNOWLEDGMENTS: please see standard acknowledgement in H.E.S.S. papers, not reproduced here 
due to lack of space. 
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